function[FSPM] = addspm(FSPM,SPM),

    nnew = length(SPM.spectra);
    nold = length(FSPM.spectra);
    k=1;
    for i = nold+1:nold+nnew,
        FSPM.spectra(i)=SPM.spectra(k);
        FSPM.spectra(i)=SPM.spectra(k);
        k=k+1;
    end;
